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Abstract 

Based on [12], we introduce a bijective correspondence be- 
tween first order differential calculi and the graph structure of 
the symmetric lattice that allows one to encode completely the 
interconnection structure of the graph in the exterior deriva- 
tive. As a result, we obtain the Grassmannian character of 
the lattice as well as the mutual commutativity between basic 
vector-fields on the tangent space. 

This in turn gives several similarities between the Clifford 
setting and the algebra of endomorphisms endowed by the 
graph structure, such as the hermitian structure of the lat- 
tice as well as the Clifford-like algebra of operators acting on 
the lattice. This naturally leads to a discrete version of Clifford 
Analysis. 

MSC 2000: Primary 39A12; Secondary 30G35 ,39A70 ,06D50 
keywords: Discrete differential forms, symmetric lattice reduction, hermi- 
tian structure. 

1 Introduction 

There are many reasons for studying discrete structures in mathematics and 
physics. 

As it was shown in a series of papers [HJ EE E], one of these reasons is 
the numerical treatment of problems related to potential theory and boundary 
values problems, where the development of discrete theoretical counterparts of a 
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continuous operator calculus leads to well-adapted numerical methods. Another 
reason is that lattices in physics exhibit some noncommutative geometric nature 
and correspond to a regularization process in field theory as well as a naive 
approximation of the topology of space or space-time, that may at high energies 
exhibit a topology different from that of the continuum. 

Although noncommutative geometry has deep roots into quantum mechan- 
ics, this notion has been introduced by A. Connes in his extension of the calculus 
of differential forms and the de Rham homology of currents [8]. 

Among the first implications was the construction of a Yang-Mills-Higgs 
theory employing the C(X) © C(X) and later (H © C) <g> C(X) algebra, where C 
and H denote the field of complex numbers and the skew-field of quaternions, 
respectively, and C(X) is the algebra of functions over X [9]. Since then, there 
exists a growing interest among theorists in studying noncommutative geometry. 
A major reason for this is that noncommutative geometry is ultimately related 
to quantum groups |27J. The latter are connected with some important aspects 
of physics, such as quantum spin chains, conformal field theories, quantum 
integrable models and so on. 

There also exists another proposal by A. Dimakis et al. jTT] in which the 
coordinates are kept commutative but there exists non-commutativity between 
the coordinates and their differentials. In this particular case the continuity is 
lost and the space acquires a canonical lattice structure. 

Regardless of the last viewpoint one looks to discrete structures as a cer- 
tain kind of differential calculi over discrete sets. First order differential calculi 
on discrete sets were found to be in bijective correspondence with graph struc- 
tures [12] . where the vertices of a graph are given by elements of the set and 
neither multiple edges nor loops are admitted. In particular, this supplies the 
elements of the set with neighborhood relations. 

The power of differential calculus rests mostly on its intrinsic character and 
on the algebraic structure which is the Grassmann algebra. In spite of their 
powerful nature, Grassmann algebras do not incorporate some concepts that 
are crucial in physics like the concept of a spinor and of a Dirac operator. This 
is one among many reasons why the correspondence between Clifford algebras 
and lattice structure should be established a priori. 

It is well known that Clifford algebras can be defined to lattice structures in 
terms of the cup algebra of simplicial homology theory. The idea was therefore 
to start from the algebra of endomorphisms of the vector space of cochains, and 
look for a natural algebraic structure in this space just like for Clifford algebras 
in the continuous case [29J. This approach, however, has some limitations, 
particularly in relation to lattice gauge theories, where the forward/backward 
differences are replaced by the symmetric differences \{d^J + d^). In 
particular, we want to write a lattice version of the Dirac operator which splits 
a lattice version of the Laplace operator. 

One way to overcome this problem was proposed by Wilson in [30] , by adding 
an extra term to the lattice version of the Dirac operator such that an extra 
fermion acquires mass and of a order of the cut-off. Another way was recently 
proposed by Faustino, Kahler, and Sommen in [17] using the splitting of the 
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standard Clifford basis vectors ej in ej = e^- + such that Dirac operators on 

lattices are constructed by using superpositions of the type e+d^ + ejd^- 1 . 

Although both approaches seem to be very promising tools, however, they 
have some limitations, particularly in relation to the lattice structure. The first 
in a certain sense requires the usage of second-order operators to define a Dirac 
operator while the second at first seems to be quite artificial. We will show 
that the second approach comes natural as a special case from the hermitian 
structure of the lattice. In our opinion, a better geometric understanding about 
the concept of Dirac operators is crucial to the formulation of well-adapted 
theoretical approaches on lattices, namely the theory of discrete monogenic 
functions. 

This paper is organized as follows: first we will introduce the discrete differ- 
ential geometry setting needed to describe the notions of exterior and interior 
product in terms of differential forms and dual connections between differential 
forms and vector fields on the lattice, respectively. Such an axiomatization can 
be found in [14] and was inspired by the works of Vaz [29], Sommen [28J and 
Brackx, Delanghe, Sommen [5]. 

Next we associate an un-oriented graph with a differential calculus on Z n by 
truncating most of the non-local links except the nearest neighboring links (i.e. 
a symmetric lattice reduction). This will be fundamental to prove some well- 
known properties of Grassmann algebras and vector-fields, namely the exterior 
product rule of differential forms and the mutual commutativity between vector 
fields. 

Afterwards we introduce the necessary endomorphisms such as the interior 
and exterior product which enables us to establish some similarities between 
the symmetric structure of the lattice and complex Clifford algebras C2n • The 
resulting Clifford structure will be compared with the Clifford-like structure on 
the symmetric lattice introduced by Kanamori and Kawamoto in [23] . 

On the last section, we introduce the Dirac operators and the corresponding 
vector variables on symmetric lattices using the hermitian setting (see, for in- 
stance, [26], [2] and [4] and references given therein). We do not claim that the 
obtained construction of discrete Dirac operators is new. It can also be found 
in the works of Vaz [29J, Kanamori, Kawamoto [23] and Forgy, Schreiber [19] , 

We also establish some basic intertwining relations between these operators 
which allow us to derive the Euler operator and also the discrete counterpart of 
homogeneous polynomials on the symmetric lattice. The resulting notion of dis- 
crete homogeneity will be compared with the hermitian homogeneity introduced 
in 0]. 

2 Discrete differential geometry 

In this section, we present a short resume of the well-known notions of discrete 
differential geometry following the works of Dimakis and Miiller-Hoissen pre- 
sented in [HI [131 H3]. The same approach can also be found in [23] as well as 
in [19]. 
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2.1 Universal Differential Algebra on a Lattice 

We consider lattice functions as maps from a lattice point I to a complex number. 
These functions jointly with pointwise addition and multiplication constitute an 
unital algebra A which is associative and abelian, and possesses unity 1 € A 

The universal differential calculus corresponds to the pair (d,A*„4), where 
A* A is a Z— graded associative algebra (over C) 

r>0 

where A r A are A— bimodules, (i.e. can be multiplied from the left and right by 
elements of ^4), and d : A r A — > A r+1 A a C— linear map which is nilpotent and 
obeys the graded Leibniz rule: 

d(dw r ) = 0, (1) 
d(Lu r uj s ) = (d.Lu r )uj s + (— l) r cj r dcj s , (2) 

where to r € A r A and lo 6 A*„4. 

From these definitions it can be easily seen that 

dl = (3) 

Assuming that A°.4 := A, we will furthermore exclusively consider differen- 
tial algebras A*„4, i.e. we require that d generates the spaces A r A. 
Let £ be a denumerable set and / a C— valued functions over C: 

/:£-C, l^f(l) = fi. 

Then the algebra is generated by the set of discrete delta- functions {b;}; 

bj(m) = Sim, 

where 5i m is the standard Kronecker symbol. 

Indeed, any function / G A can be expanded as / = J2i fl^l- This also gives 
rise to the following properties: 

b;b m = <5; m bi, 

i 

which reflects the pointwise product of functions and assures the completeness 
of the basis. From the above properties we obtain 

b;db m = -db;b m + ^ m db m (4) 

E db < = ° ( 5 ) 

I 

by using (2| and ©. 
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Using these discrete delta functions A r A can be written in terms of the 
elements 

_ J b m odb m i . . . db mr -idb m r, m° ^ to 1 , . . . , to' " 1 ^ m r 
°m°,m\...,mr-\ m r - | q, otherwise 

From ([4]) we obtain, 

t ) m D ,m 1 ,...,m r_1 ,m r = ^ m ob m l jm r-lb m r = b m odb m l . . . db m i-b m ^ 

Therefore, all the elements b m o )7n i ) ... jm r- i )Tri «- are built by concatenating el- 
ements of the form b m j m3 +i : 

^m .m 1 , . . .,m r_ 1 ,m r b m . . . b m r- 1 j77l r . (7) 

This naturally leads to the d— action 

r 

db m o _ _ _ .,m r ~ 1 ,m r ^ ^ ( 1) t) m m l _ _ . ? 77is- 1 J^ m s m r , (8) 

which represents the action of d on an r— form. 

All this shows that b^^i^^r-i^r describe simplicial r— paths which are 
related to the connectivity of the discrete space and assures the nil-potency of 
the discrete differential operator d, when acting on (J7j) 

d(db m o )m i ! ... jm r-i im r) = 

In particular, every 1-path, described in terms of non-vanishing elements 
b; m , corresponds to an edge connecting nodes I and to. This induces a graph 
structure on the set C 

In order to handle the interconnection of points of C, it is helpful to consider 
some additional structure: To make use of the additive structure of the lattice, 
it is sufficient to assume that the set C is equipped with a group addition 
I + m = p, which allows to go from I in the direction to and arrive at p on 
the discrete space and a group inverse — m defined by p + (—to) = I which 
corresponds to the direction to oriented in the opposite side. This particular 
class of graphs are Cayley graphs, i.e. discrete groups. 

The left action of a Cayley graph on the lattice will be denoted by the 
translation operator Tj defined by 

Tlf = J2fm+l^m = ^fmbt-rn- (9) 

m m 

Then every node / £ C gives rise to the 1— forms 

e l = J2*> m ,m+i, if£ m e> mb ^o, (io) 

m 

which are obviously left invariant in the following sense 

T^© — K - )• ^ ^ b rn,m-\-l — / ^ ^m-\-p,m-\-p-\-h 



(6) 
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A further interesting property of the forms 8' (and indeed of the discrete calcu- 
lus as a whole) is that they do not commute with functions, but instead induce 
translations on them: 



(11) 



b;e m , = e m b l+m 

Q l f = (T ; /)9<. 
Using this quantity we can rewrite ((7j) as 

v i r\m} —ni" at, 

Om° ,m 1 , . . .,m r ~ 1 ,m r *^mP ^ . . . \) 

In terms of 1-forms B mJ every element oj r G A r A is given by 

r 

= £ £ F m\m\...,mrQ ml - m ° ■ ■ ■ e" 1 '"™" (12) 
j=0 mi 

where 

^m°,m 1 ,...,m r /m ,m 1 ,...,m r ^m ^ *^-* (1^) 

Let us finally remark that the 1-forms (JTOj) allow to write the exterior deriva- 
tive in terms of the the suggestive action 

df = J2 M*>m,i - b,, m ) = J2( T if ~ /) 9 '- ( 14 ) 

l,m I 

2.2 Vector fields and dual connections on a lattice 

Let T denote the dual space of A 1 A as a complex vector space. Let {<9 m ,i} be 
the basis in T, dual to the basis {b/ m }, i.e. 

(bj, m , d p , q )o = Si rq 5 mtP . (15) 

T is turned into a corresponding (left or right) *4-module by introducing the 
left and right actions 

(a,f-X} Q = (af,X) , (a,X ■ f) = (fa,X) . (16) 

As a consequence we have 

b; ■ dp^n — $m,pdp,riii Qp.m ' b/ — Sp,ldp,m- 

Thus, any element X 6 T can be uniquely decomposed into 

I. Ill 

where the summation J2i m runs over an m e ^ f° r which there is an edge 
from Z to to associated to A 1 A 
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Now we introduce a duality contraction (•, •) on A 1 A as a right A— module 
and T as a left A— module by putting 

(bi, m ,X) :=bi(bi tm ,X) . (17) 

Then we have 

(fa,X-h)=f{a,X)h, (a,f-X) = (af,X). (18) 
Moreover, the elements of T become operators on A viz 

X(f) := (df,X). 
Using the Leibniz rule for d, one proves 

X(fh) = fX(h) + (h-X)(f). 

The space 

Ti:={X-hi : XeT} 

may be regarded as the tangent space at I G M, dual to the subspace 

A} A := biA 1 A = {bif : f G A 1 A} 

with respect to the duality contraction (■, •}. The space A; 1 ^ may be regarded 
as the cotangent space at I G M. 

As a consequence of this construction, {d m ^} m is a basis of % which is dual 
to the basis {b; jTn } m of A; 1 A. 

We thus have the following decompositions 

a 1 ^ = E0 A M, 

leM 

f = E0^- 

3 Hermitian and lattice structure 

3.1 Discrete differential calculi and graph structure 

In the last section we introduce universal differential algebra which is the lattice 
counterpart of a differential structure on the ordinary cases. However is some- 
times relatively far away from the concrete lattice formulations, since the the 
construction of b m j is lacking, at least initially, the neighboring structure. 

In order to obtain an appropriate lattice formulation, we need to make a 
reduction of edges by truncating most of the non-local except for the nearest 
neighboring nodes. Indeed this can be obtained from the first order differential 
calculus as a quotient A 1 A/J l , where J 1 is a submodule of A 1 ^, generated by 
non- vanishing elements of the form bm, p - 
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This kind of differential calculi appears frequently on literature. In particu- 
lar, it seems that most bicovariant differential calculi on quantum groups have 
this property, see e.g [22]. There are also several examples reduced first order 
differential calculi on discrete groups and cellular networks (see e.g. [121 [141 125] 
and the references given there). 

In our case, we will consider the type of reduction proposed by Dimakis and 
Miiller-Hoissen in [12J for the symmetric lattice. 

From now on, we will consider the case of our surface M being the Euclidean 
space R n partitioned in n— dimensional simplicial complexes, whose vertices 
form a n— dimensional lattice C, isomorphic to Z n . 

We start to consider a simplicial complex labeled by the set of vectors 
{ Vj : j e [n] U [<}, where [n] = {1, 2, . . . , n}, [n]' = {1', 2', . . . , n'}, ' : 
[n] U [n]' — > [n] U [n]' is an involutive permutation which maps j ^ j 1 and j' i— ► j 
and j h-> Vj , j' i ^ -Vj 

This lead to the lattice reduction 

This reduction implies that there are 2n nonzero differentials 0- with 

6 Vj = bm,m+v 3 -, 3 G [n] U [n]'. (20) 

Let us remark that, by construction, the vectors and Vf have equal lengths. 
Then the translations m + Vj and m + v^v are symmetric to each other with 
respect to the hyperplane TL which contains the point m of the coordinate system 
in E n and which is perpendicular to the edge which links the nodes m + Vj and 
m + Vj> . 

This means that in our lattice reduction, the edges between the nodes m, m+ 
Vj and m, m + Vj> are being kept and thus links are symmetric with respect to 
their orientation. 

Therefore, the important automorphisms on A* A leaving the symmetric lat- 
tice reduction invariant are the main involution defined by 

(uj r uj r y = u' r uj' s , 

;/'v v (-> v . . . b v • i' = F Ju ... Jr (e v ny . . . (9 V -)', (21) 

the reversion given by 

(LU r LU r ) — UJ s UJ r , 

v ■ ..« v ; = ;B V * • (e v «)~F Vii> ..., v , r , (22) 

(e v o =-e v 3 ', 

and, finally, the f — conjugation defined by 

(<jj r u) r y = uj\ uj,, 

i/ v v B v • . -O^y = (9^)t . . . (Q^)W Vjii ... tV . r , (23) 

(e v o f = -e v ^', 
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(we remark that F Vji> ... tVjr = f Vj jVj bm stands for the complex conjugation 
on A). 

As explained in [12], first order differential calculi on a finite set are in bijec- 
tive correspondence with digraphs with at most a pair of edges connecting two 
vertices. Such a graph is characterized by its adjacency matrix which is a square 
matrix Q such Gm, P = 1 if there is a edge from to to p and Gm,p = otherwise. 
Therefore, we can associate with each non-vanishing bm lP of some differential 
calculus (A* A, d) an (undirected) edge connecting the nodes (vertices) m and 
p and hence the adjacency matrix Q is assigned by the sum of all 1-forms 

G=E E b m,m+v, = E 6VJ - ( 24 ) 

m££ je[n]U[n]' je[n.]U[n]' 

The first order differential calculus then corresponds to a complete graph 
where all vertices (the elements of C) are connected to each other by a pair 
of edges, whose interconnection structure is completely encoded in the exterior 
derivative d. Indeed, from properties and Ijlip . it follows that our first order 
differential calculus is inner since G acts like d when commuted with functions 
/ e A, i.e. 

[G, /] := J^/mCGbja-bjaG) 

= E E /m(® Vj bm- b Za 0V ') 

me£ je[n]U[n]' 

^ ^ ^ ^ /m (bm-vj ,m t*;m.m+Vj ) 

me£ je[n]U[n]' 

= df (25) 

To the nearest-neighbour-nodes being encoded in the 1-form b^m+v,-! we 
associate a tangent vector field dm+ Vj t! n G %n ■ 

From the d— action lfl4|) and the contraction constraints lfl5|) and lfl6|l . we 
can consider dm+vj ,m £ as a partial difference action at node to: 

^m+Vj,m ,/m — fm+vi /m ; foi 3 — 1, . . . , 71. 

Moreover, a basis {d Vj }je[ n ]u[n]' f° r the vector space T satisfying the duality 
relation 

(Q v *,d v *) = 5 kj (26) 
is uniquely determined by the relation 

d Vl f = ]T ^M+Vi ,m/m)b„ = T Vi / - /. (27) 
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Therefore, 

d/ = [G,/]= (^ Vj /)e Vj (28) 

je[n]u[n]' 

From the above characterization, we can now prove some results which relate 
d with the graph structure assigned by G. 

According to definition (fSJ) and property (Q]), the set of r— forms is completely 
determined by the set of edges of the graph. Indeed, the sum of all r— forms 
then corresponds to the r— power of G, 

G r = Yl e Vfc i ...e Vfc -, (29) 

{ki,...,k r }c[n]U[n]' 

which is nothing else than the r— power of the adjacency matrix Q. 

Let us take a close look at the sum of all 2— forms. On the right-hand side of 
l(29|). the summation over all indices j, k £ [n] U [n]' involves terms of the form 

_ (0v 3 @v fc _|_ QVfcQVj-j _ _ ^2 (b BB+Vjj!! , +V]+Vt + b M+Vta+Vj+Vl ) (30) 

Because of bm,m+vj+v k = the right-hand side of l(30|) in terms of the action 
of d then corresponds to 

1 (e v 36 v fc + eVkQVj) = lj2 db m , m+Vj+Vfc =0, j+ k'. (31) 

Contrary to the oriented lattice reduction, in the symmetric lattice reduction 
the left-hand side of 1|30|1 does not vanish in general (see e.g. [23]). However 
from ((3TJ) and 

n 

o = 53 db™,™ = - 53 (e v *e v i' + e v i'e v o 

the square of G, G 2 , vanishes. 

Therefore, to ensure that the term l|30p vanishes for all j, k 6 [n] U [n]' , our 
graphs should satisfy the following two conditions (see the paper by E. Forgy/U. 
Schreiber p]): 

1. G has no intermediate edges: hm, P 7^ => hm,i P = 0, 

2. G has no opposite edges: bm iPil = 0. 

Notice that on graphs without intermediate edges all the 2-paths of edges 
assigned by b Bi ; iP that connect the same two points has to vanish. So our intro- 
duced lattice reduction (fl9l) is a particular case of a graph without intermediate 
edges. 

On the other hand, on graphs without opposite edges, neither multiple edges 
neither loops are admitted. 

Therefore the above conditions leads the following theorem: 
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Theorem 3.1 For graphs without opposite edges, the exterior product rule 

{e v %e Vfc } = o, 

holds for all j, k e [n] U [n]' . 

Here and elsewhere, {a, b} — ab + ba denotes the anti-commutator bracket of a 
and b. 

By applying the coordinate definition ([28| . Theorem (|3.ip and the nil- 
potency of d ([J), we arrive at the following corollary: 

Corollary 3.1 For graphs without opposite edges, the basis elements d Vj ,j S 
[n] U [n]' , mutually commute when acting on functions A, that is 

g Vj ( d v k fj = g Vk (ff*if^ f or aU j > fc e [ n ] u [n]' and f G A. 

Let us remark that, for the product of two functions f,g € A, the difference 
action 9 VJ satisfies the product rule 

d^{fg) = (d v if)(T Vj g) + f(d^g). (32) 

This means that the finite difference actions d Vj do not obey the ordinary 
Leibniz rule. In fact, the application of d Vj to, say, higher powers becomes 
increasingly cumbersome. 

Due to the discreteness of the formalism and, as a consequence, the inevitable 
bi-locality of d Vj there is no chance to get something as a 'true' Leibniz rule on 
this level. Nevertheless, there is a certain systematic in it, namely the product 
rule. 

We will come back to the non-Leibniz character of d Vj , when establishing 
the duality between differential forms and vector-fields. Before however doing 
that we will need to further clarify the role of d on r— forms. 

Under the conditions of Theorem 13.11 we get 

fce[n]U[n]' 

Then we have 

Applying induction on r > and the graded Leibniz rule 

d(u>iujr) = (du>i)u) r — £Ji(duv) 

together with the action of d on r-forms ([6]), lead to the following theorem. 

Theorem 3.2 For graphs without opposite edges, the actions of G and d on 
r— forms are related as follows: 

duj r = Gw r — (— l) r u; r G 
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3.2 Differential forms representations of End(A*A) on the 
symmetric lattice 

In order to introduce directions in lattices with mesh- width h > 0, we introduce 
the coordinate functions 



2J h rUjhrn 



(33) 



and split up G into G = G + G, where 

g = e;=i© v ^ g = e;=i© vj - 

Hence, by applying d on Xj, we can find the coordinate differentials dxj on the 
symmetric lattice as 



dxj = [G,Xj] = —dx- + dxf . 



(34) 



where dx^ corresponds to 



hO v ^ dx^ 



-G,Xj 



According to I133II we can express any function / 6 A as f(x) and by direct 
application of the coordinate differentials (|35| , the non-commutativity of func- 
tions and 1— forms shown in iflTj) can now be given in terms of x — (xi, . . . , x n ) 
as 



dxff{x) = f(x±hv j )d 



(35) 



To represent any r-form u> r € A* A in terms of the coordinate differentials 
dx~~, we proceed as follows: 

For two ordered subsets J' C [n]' and K C [n] given by 

J'={j[,...,f s }, K={k 1 ,...,k t }c[n], 

we put 

dxg = 1, 
dxj = dxj ± . . . dxj s 
dx+ =dx+ i ...dx+. 



Therefore, the r— form defined by l|T2|) then corresponds to 
u r {x) = ^2 Fj> K (x)dxjdx~^, 

\J'\ + \K\=r 

where Fj'k denotes the function (jT3j) indexed by J' UK and J2\j'\+\K\=r denotes 
a sum restricted to ordered subsets J' C [n]' and K C [n]. 
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For the vector fields £ T defined as the forward/backward difference 
actions 

{dt J f){x) = \(d^f)(x) 1 (d^f)(x) = -l(d^-'f){x), Vf£A (36) 

we can split the d— action on oj r as 

cLlv(x) = d + Lu r (x) — d-U r (x), (37) 
where d±ui r (x) are given by 

n 

d±u> r (x) = Yl J2 (d± j Fj, K )(x) dxfdxjdx+. (38) 

3=1 \J'\ + \K\=r 

The above framework suggests the following bi-graded algebra (i.e. a bi- 
complex) decomposition 

n 

A* A = © AM -4 

p,q=0 

where A := A 0: °.A and the discrete exterior differential maps 

d : A™ A -> A p+1 ' q A, d+ : A p ' q A -> A*** 1 A 

are defined by formula Under the conditions of Theorem GTTJ we have that 
the coordinate differentials dxf satisfy the anti-commutation relations 

{dxf,dx^} =0 Vj,fc=i,..., n , 

{<l.'-;.clr, } v,- fc=1 ,...,„. 1 ' 

and, furthermore, 

d±(d±uj r (x)) = = d + (d_w r (x)) + d_(d + w r (x)). 

The basic endomorphisms acting on A* A in an exterior way are the linear 
operators 7 ±J £ End(A*_4) defined as 

f-j : A™ A -► A^'M, w(s) i-» da?7w(as), 

7 +3 . ap.M -> A^+U, w(x) h-> dx+w(x). ^ ' 

Having defined the left and right exterior product representations, it raises the 
question how to define left and right representations for the interior products 
# i £ End(A*„4) in terms of the duality contraction (■,•). To be consistent with 
the nature of the interior product, we impose the duality conditions ^^(dxt) = 
and d ±j (dx^) = 6 jk . 

Let us proceed as follows: first since the Leibniz rule 10 is also valid for the 
exterior derivatives d± we observe 



d±{xkUJ r {x)) — dx k LJ r (x) + Xkd±Lu r (x) 
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and, hence, direct application of the duality contraction property (|26|) and the 
non-commutativity (|35| lead to 



dx*u(x),d± J 



S jk (T^u)(x)-dx 




)A ±J ')- 



(41) 



Hereby, T h 3 denotes the shift action (T h 3 uj){x) — u(x ± hVj) on the exterior 
algebra A* A. 

This suggests that the interior product operators £ End(A*„4) should 
be defined as contraction operators with shifting role opposite to the differential 
form 



There are some formulae that follow from equations <[39|l . (|40|l and (|42ll will be 
of interest, namely, the anti-commutation relations between the interior and the 
exterior product representations 7 J , 1? J 6 End(A*A) 



7 ±i (7 ± M£)) +l ±k (l ±j u(x)) = 0, 

j +j {~f- k Lu{x)) +7 (i +j lu(x)) = 0, 

^(■d ±k uj(x)) + $ ±k (ti ±j uj(x)) = 0, 

+ '(0-*w(ss)) +tf- k ({>+3Lj(x)) = 0, 

7 +3 (r i: a)(i)) +?9- fe (7+^y) = 0, 

7 ±i(0±fc w (z))+0±*( 7 ±J w (£)) = ^few(x) 



The above identities then correspond to the graded algebra like the one of 
Fermionic creation and annihilation operators |24j . 

It's now interesting to compare our approach with the approach proposed 
by E. Forgy/U. Schreiber in [l9j, where they construct a discrete differential 
calculi on causal graph complexes. 

On their terminology (see [19], page 4), the algebra of endomorphisms End (A* .4) 
would be the inner product space H(A, (*|-)). Indeed a representation of the ba- 
sic endomorphisms 1? J in terms of the inner product (-|-) already tacitly exists 
(see [19], section 3.6). However, bear in mind that these operators can be in- 
troduced independently of the inner product (-|-) and the behavior is invariant 
under diffeomorphisms (c.f. PQ). 

Moreover, our framework gives a quite natural way to describe the alge- 
bra End(A*„4) as a canonical realization of the Fermi algebra in a metrically 
independent, via duality arguments. 

We are now in conditions to explore the correspondence between differential 
forms and Clifford algebras on the symmetric lattice. 



: AP<*A -> hP-^A, uj{x) 1 — y 
&+i : A p q A -» AP'i-i-A, uj(x) i-> 



{Ti^){x),d~ h 



(42) 
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3.3 Clifford algebras and hermitian structure of the sym- 
metric lattice 



As it was shown in [3, 20, 28J Clifford algebras can be defined in several different 
ways. One of these ways is as a subalgebra of the algebra of endomorphisms of 
the exterior algebra. 

In order to get this correspondence for the symmetric lattice, we proceed as 
follows: 

First of all, notice that from the commutation relations (|43|) the endomor- 
phisms : j = 1, . . • ,n, € End(A*.A) defined by 

£±j _ 7 ±i + (44) 

satisfy the graded fermionic identities when acting on AM : 

^(rM£)) + r fe (£ + M£)) = torallweAA 

Furthermore, the elements 

satisfy the graded orthogonal identities when acting on A*„4 : 

r^(r^(x)) + r^(r^u,(x)) = ±26 jk u,(x), 

T+J'(T- fe cj(x))+T- fe (T+J'cj(x)) = ° a fc A 

This clearly suggests that T ±J behave like the generators of the real Clifford 
algebra of signature (n, n), 

For the symmetric lattice we get the following interesting features of the 
coordinate differentials dxj — Y _J (1) and dTj := dx+ + dxj = T +J (l): 



(dXj)' = 


-dxj, 


(drj)' 


= dTj 


(dxj) 


= dxj, 




= 


(dxtf = 


-dxj, 


(dT,) f 


= dTj. 



Hence dxj behaves as a real while dr, behaves as pure imaginary with respect to 
the involution, reversion and f— conjugation, respectively. It then turns out that 
the automorphisms lf2"Tj) - ([23"l) play the same role as the automorphisms on the 
complex Clifford algebra C2« according to the symmetric nature of the lattice. 

Since the real Clifford algebra M„ iTl is contained in the complex Clifford 
algebra Ca n = C(££)Ro,2n as a special subalgebra [3], the isomorphism between 
the End(A*_4) and £-2n = C ® Ro,2n is thus obtained through the identification 
ej «-> T - -' and ej+ n «-> iT +J . 

The corresponding Witt basis for C2« is given by 

fj = — (ej — ie n+ j), fj = — -^(e j + ie n+ j) 
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and satisfies the following anti-commuting identities 

Grassmann identities: {fj, %} = = {ft, f£}, 
duality identities: {fjjfD = $jk- 

Moreover, we can identify it with the set of endomorphisms 

■ J = l,-..,n} cEnd(AM) 

viz f , <-> Z+i and ft <-> 

All the above identifications clearly suggests the hermitian Clifford setting 
as the natural multi-vector setting to develop a discrete framework on the sym- 
metric lattice. 

It is now interesting to compare our approach with the approach proposed 
by Kanamori and Kawamoto in |23j. 

We would like to remark that the initial construction proposed in subsection 
13.11 is analogous to the approach proposed in that paper. However, while the 
authors define a Clifford-like algebra on the symmetric lattice by introducing 
a new Clifford product in a rather "artificial" way, we avoid this problem by 
passing to the algebra of endomorphisms. 

While in the approach of Kanamori and Kawamoto in [23J there is no restric- 
tion to obtain Clifford products which are associative and distributive, however, 
that construction is metric-dependent, since the Clifford product can be indefi- 
nite. 

In our approach, we note that a priori no problem in relation to the question 
of associativity and distributivity on the lattice will appear since the algebra of 
endomorphisms of a given space equipped with the standard sum and a product 
defined by composition is obviously associative and distributive. On the other 
hand, the structure of the algebra of endomorphisms lead also to a construction 
of Clifford products which are metric-independent in the sense that they only 
depend on the duality between the tangent and the cotangent space. 



4 Dirac operators and vector variables on the 
symmetric lattice 

The main objective in this section is to show some similarities between the 
hermitian setting and the symmetric structure of the lattice. 

Our starting point again is the definition of d± and d = d + d . 

From Corollary 13. II we know that all forward and backward differences mu- 
tually commute when acting on A* A 

Furthermore, they are interrelated by the translations (T h ^uj){x) = ui(x± fiVj) 



16 



Using f35|) and <f38|l . we can thus write the exterior differentials d± in the 
form 

n 

d ± =£dz±df. 

3=1 

Let us now introduce the symmetric and skew-symmetric difference operators 
and V^, respectively, as: 

Then the operator d defined in l|37p corresponds in terms of the coordinate 
differentials dxj and drj to 

d = d+ d 

n 

3=1 

n 

= Y^^-'M ■ 'dr,Vi. 

3=1 

The above identities suggest the introduction of the following operators acting 
on End{k*A) 

n 

d± = (46) 

3=1 
n 

9 = 'Vi • iT-^l, (47) 

3=1 

where £ J and T ±J are the basic endomorphisms defined in f44|) and (|45|) . re- 
spectively. All the geometry of the symmetric lattice is now encoded in the 
operators d± and d, which are nothing else than the hermitian and the orthog- 
onal Dirac-Kahler operators acting on symmetric lattices, respectively. 

The operator defined in lj47|) is closely related the Dirac operator intro- 
duced by Kanamori, Kanamoto in [23] (see formula (4.10), page 21) and Forgy, 
Schreiber in [19] (see formula (5.34), page 71). 

Here we are interested to establish a correspondence with the Hermitian 
setting proposed in [I], our departure point will be the operators (06]) instead of 
the operator (|47ll . This is indeed the main difference between our approach and 
the approaches proposed by Kanamori, Kanamoto in [23] and Forgy, Schreiber 
in [19]. 

From the correspondence fj <-> £ +J and fj «-» we see that d± is the 
lattice counterpart to the hermitian Dirac operator and its conjugate on C.2 n 

d - V" f-d~ j 
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where d^ 3 — | ^V^ =F *V^J are the discrete counterparts of the classical Cauchy- 
Riemann operators and their conjugates. Moreover, from the correspondence 
ej <-» Y - -', ej +n <-* iT" 1 "- 7 we see that 9 is the lattice counterpart of the discrete 
Dirac operator on C.2n 

n 

dx = 5>,V£+e,+„Vi (49) 

i=i 

Similarly, the lattice counterpart —i(d + + <9_) of the discrete Dirac operator on 
C 2n is given by 

n 

d x \ ■■= ~i(d z + 3\) = e M - e ^+" v i (5°) 

There are some formulae holding for the operators (J48J) , (J49J) and f50|) that will 
be of interest in the hermitian setting, namely 

1. Isotropy condition: &l — = (<9|) 2 

2. Orthogonality condition: {<9x , <9x|} = 

3. Star Laplacian splitting: 

• Using Dirac operators: d\ = — X)j=i ®h J ®h = 

• Using hermitian Dirac operators: {<9 Z , 9]} = X)j=i d^ 3 d^ 3 

Having established the correspondence between hermitian Dirac and Dirac 
operators on the symmetric lattice and discrete Dirac operators on the Clifford 
algebra, we arrive now at the question how to make the correspondence between 
coordinate vector functions on the symmetric lattice and vector variables on the 
Clifford algebra. 

Let us remark that the finite difference action d^ 3 acting on A* A satisfies 
the product rule 

8* (x k (T* k u;)(x)) = 6 jk T* k (T± 3 u>)(x)+x k T* k (d± 3 u)(x) 

= 6 jk u(x)+x k T* k (d± j cj)(x). (51) 

which establishes the duality between the finite difference operators d^ 3 and 
the "formal" coordinate functions XjT^ 3 on End(A*.A). 

We also note that the coordinate variables xfT^ 3 mutually commute, when 
acting on functions on A*^4 

*jTit j {*kT£ k u(x)) = x k T± k { Xj T± 3 u>{x)). (52) 

Furthermore, the commutative relations (f46|l . ([52| together with the duality rela- 
tion lj5T]l endow an algebraic representation of the Weyl-Heisenberg algebra [24] , 
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where the "formal" coordinate functions XjT h 3 represent "creation" operators 
dual to the "annihilation" operators d h ° . 

From the previous relations, (f43~)) . and definition l(42|). we immediately ob- 
tain the following commutation relations between d^ J and 7 ±J when acting on 
End (AM) 

dt j {^ k {u{x)))=^ k {dt j {u{x))). 

Moreover, Clifford-like operators on the symmetric lattice are therefore encoded 
in the following algebra of endomorphisms 

Alg{d^,d+ j ,x j T+ j ,x j T^,t j ,r j -J 1 "}• 

where £ ±3 6 End (AM) defined in lj44j) satisfy the fermionic relations (J45j) . 

Let us define formally Zj '■= XjT^ J , zj := XjT^ J as the complex variables 
Zj and their conjugates Tsj, respectively. Defining the hermitian vector variable 
z = Yjj=i fj z j an ^ its hermitian conjugate = J2]=i fj^J> the Clifford vector 
variable X associated to dx takes the form X = z — z* while the Clifford vector 
variable X\ associated to d X \ takes the form X\ = —i(z + z^). 

Hence, the following "formal" vector variable identifications naturally follows: 

z — E"=i t +j {*iT+% zt ^ E « =i r j {x . V)) 

^ Z^j=l 2 T 4 - 1 2* ' 

There are some formulae holding for the vector variables X,X\, z, and that 
will be of interest in the hermitian setting, namely 

1. The isotropy condition: z 2 = = (z^) 2 

2. The orthogonality condition: {X, X\} = 

3. Square variable splitting: 

• Using Clifford vector variables: 

n 

X 2 = -J2(xjT^)( Xj T^) = (X\) 2 

• Using hermitian vector variables: 

n 

{z,zt } = ^ (a;j T^)(x,V) 
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By Clifford geometric product OH], 

ab = a»b + aAb 

one can also introduce the discrete hermitian Euler operators as E z — 2z • 
d z , E z -\ = 2z^ • d\ and the hermitian Gamma operators as T z — z A d z , T z t — 

z f A dl 

The Weyl-Heisenberg character of the operators and XjT^ then lead to 
the following elegant formulae similar to the intertwining relations concerning 
the operator d z ,d\,z,z\E z ,E z ^,T z and T z \ in [2111]. 

{z,8 z } =j3 + E z , [z, 8 Z ] = -13 + T z 

{zi,di} = (n-0)+E zU [«t,fl]]=-(n-/3)+r, t (53) 
{z\d z } = Q, {*,0]} = o 

where = Y^j=i Vjh denotes the spin-Euler operator. 

Analogously, the link between the hermitian operators E z ,E z i , T z , T z t and 
the traditional Euler and Gamma operators Ex = —X • dx,E x \ — —X\ • 
dx\, Tx = —X A dx and r^i = — X\ A d x \ can also be obtained: 

E X =E Z +E Z , =E X] , 

t x = r 2 + r, t - 2(zt A d z + z A d% 

T x = T z + r, t + 2(zt A d z + z A fit). 

Let us take now a close look at the concept of discrete homogeneity on C2 n • 

As homogeneous polynomials are expected to be C2 n — valued eigenfunctions 
of the Euler operator corresponding to the eigenvalue given by the degree of 
the polynomial, in the hermitian setting it still makes sense to define discrete 
homogeneous polynomials of degree (p, q) as solutions of the coupled eigenvalue 
problem 

E z [Rp^{z, 2+)] = pR Pt q(z, zt) 

(54) 

E z \ [Rp, q (z, = qRp, q (z, zJ) 

and discrete hermitian monogenic homogeneous polynomials of degree (p, q) as 
the solutions of 1(54]) satisfying the zero hermitian Dirac constraints: 

d z R p . q (z 7 z^) = = d z tRp. q (z,z r ) (55) 

A further consequence of l|53"ll - lf55|) . is that any discrete hermitian monogenic 
homogeneous polynomial of degree (p, q) is an eigenfunction of the hermitian 
Gamma operators, namely Y z [R p _ q (z, z^)] = —pR p , q (z, z*) and T z t [R p , q (z, z')] = 
-qR Pj g(z, z^), 

It is interesting to see the similarities between the solutions of the coupled 
eigenvalue problem 1(54]) and the concept of discrete homogeneous polynomials 
introduced in our previous paper [16J. 
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In the terminology of that paper, = Xw=i x j^h^ m ^ ne forward/backward 

Euler operator while (x)± are the multi- index factorial powers satisfying the 
eigenvalue property 

£±(# = H(#. 

On the other hand, discrete hermitian Euler operators coincide with the for- 
ward/backward difference Euler operators, namely E z — E^ and E z \ = E^ . 

From the above relations, a solution of 1|54|1 is explicitly given by linear 
combinations (possibly C2 n — valued) using monomials of the type 



with \a + 1 = p and \oT \ = q. 

The crucial difference between the homogeneous polynomials formulated 
in [4] and the solutions of the coupled problem f54|) is that they do not sat- 
isfy 

However, they form a Sheffer sequence of polynomials with respect to the opera- 
tors and XjT^ , since (x)^' are basic monomials (i.e. d^Xk — §jk> (x}±^ = 
1 and (0)^ = 0.) and satisfy the monomial principle 



A further consequence of the above relations is the Rodrigues formula 

(x)i Q) = { Xl T^r{ X2 T± 2 r ■ ■ ■ (x n T± n r-i. 

Moreover, (z p 1) ((zt)? 1) is a discrete homogeneous polynomial of degree (p, q) . 

The above formulae together with the Weyl-Heisenberg character of d h ^ and 
XjT^ 3 clearly suggests a correspondence between the R— polynomial algebra 
generated by (x)±^ and the Bose algebra. In fact an isomorphism between both 
algebras already tacitly exists [6]. 

This is the starting point for constructing discrete versions of hermitian 
Fischer decompositions in terms of discrete hermitian monogenic homogeneous 
polynomials, i.e. polynomial solutions satisfying ([54)1 and (|55| and, moreover, 
generate Hermite polynomials as an Appel sequence associated with the orthog- 
onally shift-invariant Weierstrass operator [6j. This will be one of the main 
topics to be studied on the forthcoming paper p.8] from the umbral calculus 
point of view. 
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